Upper Limit on Gravitational Wave Backgrounds at 0.2 Hz with Torsion-bar Antenna 
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We present the first upper limit on gravitational wave (GW) backgrounds at an unexplored 
frequency of 0.2 Hz using a torsion-bar antenna (TOBA). A TOBA was proposed to search for 
low-frequency GWs. We have developed a small-scaled TOBA and successfully found f2gw(,f) < 
4.3 X 10^^ at 0.2 Hz as demonstration of the TOBA's capabilities, where r2gw(/) is the GW energy 
density per logarithmic frequency interval in units of the closure density. Our result is the first 
nonintegrated limit to bridge the gap between the LIGO band (around 100 Hz) and the Gassini 
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Introduction. - 

The standard cosmology predicts not only the cosmic mi- 
crowave background (CMB) from the last scattering sur- 
face, but also gravitational wave backgrounds (GWBs) 
as the result of processes that take place very shortly 
after the big bang [if. Astrophysical GWBs could also 
emerge from the superposition of a large number of un- 
resolved sources 0. The proper frequency characteristic 
of the cosmological GWBs is determined by their gen- 
eration mechanisms and the state of the Universe when 
the wavelength of the GWBs crossed the Hubble horizon. 
The frequency characteristic of the astrophysical GWBs 
depends on the motion scale of GW sources. Measure- 
ments of GWBs in various frequency bands enable us to 
separate the GWBs according to their origins, and then 
reveal how the Universe evolved from its very early epoch. 
Therefore, GWB detection and characterization is one of 
the greatest challenges in not only GW experiments, but 
also cosmology and astronomy. 

A number of experiments have been performed to de- 
tect GWBs or constrain rigw(/) at several frequencies, 
where Qg^{f ) is the cumulative energy density of GWBs 
per unit logarithmic frequency, divided by the critical 
energy density to close the Universe. At around 100 
Hz, the Laser Interferometer Gravitational Wave Obser- 
vatory (LIGO) has found f^gw(/) < 6.9 x 10"^ at the 
95% confidence level [3] . A pair of synchronous recycling 
interferometers has set fJgw(/) < 1.2 x 10^^ for 100 MHz 
GWBs |4|. At 907 Hz, a cross-correlation measurement 
between the Explorer and Nautilus cryogenic bar detec- 
tors has placed rjgw(/) < 120 At 10"^ _ iqs Hz, 
the Gassini spacecraft has established an upper limit by 
using spacecraft Doppler tracking [il . Based on the fluc- 
tuations in the pulse arrival times from PSR B1855-I-09, 
an upper limit has been found in the frequency band 
10"^ - lO"'^ Hz [l|. Measurement of the CMB at large 



angular scales indicates an upper limit at very low fre- 
quencies (10"^® - 10"^^ Hz) 0. In addition, the inte- 
grated ilgw = / Jlgw(/)d(ln/) is indirectly constrained 
by the helium-4 abundance resulting from big-bang nu- 
cleosynthesis (BBN) [l[ and measurements of the CMB 
and matter power spectra [H. 

Currently, the low-frequency band at the range of 
0.01 — 1 Hz is still unexplored, although several low- 
frequency GW antennas are being proposed and devel- 
oped, -e.g., the Laser Interferometer Space Antenna Q, 
DEGi-hertz Interferometer Gravitational wave Observa- 
tory [10 1, and Atomic Gravitational wave Interferometric 
Sensor llj . 

In the previous Letter, we have proposed a torsion- 
bar antenna (TOBA) to search for GWs at 0.01 - 1 Hz 
[l^ . For example, it is realistic to achieve a GW strain- 
equivalent noise level of ~ 10-18 Hz-i/2 at 0.1 Hz, even 
with a ground-based configuration. In this configuration, 
the observable range reaches about 10 Gpc for 10^ Msun 
black hole events. To evaluate TOBA's capabilities for 
low-frequency GW observations, we developed a small- 
scaled TOBA and performed a short observational run. 
This Letter presents the first upper limit on GWBs at 
0.2 Hz obtained from the small-scaled TOBA as demon- 
stration of the capabilities, and discusses the future strat- 
egy of a TOBA. 

TOBA.- 

A TOBA differs from conventional ground-based GW an- 
tennas in its fundamental sensitivity to GWs below 1 Hz. 
A TOBA consists of two rotational sensors and two bar- 
shaped test masses arranged parallel to the x-y plane and 
orthogonal to each other (see Fig. [1]). Each test mass is 
suspended from its center; therefore, it behaves as a free 
mass in the rotational degree of freedom around the z 
axis. When x -polarized GWs of /i^ (t) (/ii2(i) — h2i{t) — 
hy,{t) and h^j = ^ (1, 2), (2, 1))) pass through a 
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FIG. 1: Principle of a TOBA. Two orthogonal test masses 
feel the tidal forces by incident GWs. 



TOBA, tidal forces by the GWs will appear as angular 
fluctuations of the two test masses. 

The angular fluctuation 6* of a test mass from its orig- 
inal position obeys the equation of motion fl3] : 



(1) 



where /, qi2 and /res = 1/ (27r) ■^/T/k are the moment of 
inertia of the test mass, the dynamical quadrupole mo- 
ment {qi2 — (721 and qn — —922, see [ISl]), and the rota- 
tional resonance frequency, respectively. Here, 7 and n 
are the damping constant and the spring constant around 
the z axis, respectively. A tilde (~) denotes the Fourier 
transformation. 

The equation of the motion of another test mass is also 
written as 9' = — (7i2/(2/)ft,x (/), where 9' is the angular 
fluctuation of the test mass. The differential fluctuation 
A9{t) (= 9{t) - 9'{t)) is expressed as. 
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GWs can be detected from 9 with the single-mass config- 
uration or A9 with the differential-measurement configu- 
ration at /res ^ / ^ //(27r)7. In general, the resonance 
frequency /res in a torsion pendulum is as low as a few 
mHz [14|. Thus, a TOBA can have a fundamental sensi- 
tivity in low frequency band (0.01 - 1 Hz). 

Small-Scaled TOBA- 
We developed a small-scaled TOBA as shown in Fig. [21 
The small-scaled TOBA mainly consists of a single test 
mass and a interferometric rotational sensor. As the first 
step, we applied the single-mass configuration, though 
we had to abandon the common mode noise rejection ob- 
tained from the differential measurement. Other unique 
equipment is a magnetic suspension based on the pin- 
ning effect of a type-II superconductor. This suspension 
can potentially provide large suspension forces with the 
low spring constant n and the low damping constant 7 
without the drawbacks of contact 
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Our test mass has an inverted T-shape with a hori- 
zontal length of 22.5 cm, a vertical length of 19 cm, and 
a mass of 131 g. On the top of the test mass, a sym- 
metric neodymium (Nd) magnet (022 mm, tlO mm) is 



attached. The superconductor, which is made of a Gd- 
Ba-Cu-0 compound and has a critical temperature of 
92 K, is placed above the Nd magnet. When the super- 
conductor is cooled to about 65 K by a low-vibration 
pulse-tube cryocooler, the test mass is suspended by 
the magnetic forces between the Nd magnet and the 
superconductor. The shape of the test mass was cho- 
sen to bring the Nd magnet close to the superconduc- 
tor and elongate the arms as much as possible. The 
spring constant and the damping constant of our sus- 
pension system are k = 3.6 ± 2.1 x 10~^ Nm/rad and 
7 = 1.2 ± 0.7 X 10^* Nms/rad (a rotational resonance 
frequency of about 5 mHz and a quality factor of about 
1000), respectively. These values are similar to those of 
a typical tungsten-based torsion pendulum. The current 
damping constant is limited by the random collision of 
residual gas under 10"'^ Pa condition. 

We implement a laser Michelson interferometer using 
two mirrors attached on each end of the test mass. The 
angular fluctuation 9 then is measured from the optical 
pass difference Al in the interferometer. The fluctuation 
9 is fed back to the test mass using coil-magnet actuators 
for the operation of the interferometer in its linear range. 
In this configuration, GW signals are read out from the 
feedback signal, just like the conventional antennas [16 1. 

As a laser source, a Nd:YAG laser is used with a wave- 
length of 1064 nm and an output power of 40 mW. To 
compose the actuators with external coils, the test mass 
houses two samarium-cobalt (SmCo) magnets {4>1 mm, 
to. 5 mm) at each end of the horizontal arm. The test 
mass and most of the readout interferometer are located 
in a vacuum chamber, magnetic and thermal shields to 
avoid the effects of external disturbances. 

We find a GW strain-equivalent noise level h{f) of 
2 X 10^3 Hz- around 0.2 Hz (Fig. [3|). Current noise 
level is limited by environmental disturbances: seismic 
noise above 0.2 Hz and magnetic noise below 0.2 Hz. 
The seismic noise is from the unexpected coupling be- 
tween the translation of the test mass induced by seismic 
motion and Al in the interferometer. The magnetic noise 
is introduced by the coupling between the residual non- 
axial symmetry of the Nd magnet and the fluctuation of 
external magnetic fields. Little differences of the struc- 
ture between the measured noise level and the estimated 
environmental noises shown in Fig. [3l are due to measure- 
ment conditions. The peaks at 0.3 Hz, 1.2 Hz, 3.5 Hz and 
7 Hz are identified as the effect of microseismic distur- 
bances, the resonance of the rigid-body-pendulum of the 
test mass, the resonance of the platform of the chamber 
and the resonance of the lab's floor, respectively. Com- 
pression of the cryocooler also induces peaks at 3.9 Hz 
and 7.8 Hz. 

More details of instrumental status, including noise 
analysis and performance of the magnetic suspension, are 
published as individual papers 1714191. 

Upper limit on GWBs- 
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FIG. 2: Conceptual design of small-scaled TOBA. 
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FIG. 3: Measured GW strain-equivalent noise level. Es- 
timated environmental noises (seismic noise and magnetic 
noise), technical noises (sensor noise and electronics noise) 
and thermal noise are also described. 



On August 15, 2009, we performed a one-night observa- 
tional run using the small-scaled TOBA. The total length 
of the data was about 7.5 hr. 

Using approximately Toff — 3.5 hr of stable low-noise 
data, we place an upper limit on J7gw(/o) at /q = 0.2 
Hz with a bandwidth of 10 mHz. Such arbitrary data 
selection is acceptable, since we assume stationary GWBs 
in this Letter. 

The data analysis pipeline consists of estimating a GW 
energy-equivalent spectral density rioq(/o) and determin- 
ing the upper limit r2g^(/o). Supposing that GWBs are 
isotropic and unpolarized, the ricq(/o) can be related to 



f^cq(/o) 



10^ 



(3) 



Here, Hq is the Hubble constant, and its value is 70.4^}'4 
km/s/Mpc [13 . 

To obtain the estimator of r2oq(/o), the data are di- 
vided into ts = 102 A sec segments. The length of the seg- 
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FIG. 4: Distribution of ncq(/o). Green boxes are rejected 
segments. Red curve is the result of the exponential fitting. 
Blue bars are Poisson errors of the distribution. 



ment is selected such that the bandwidth is much smaller 
than the target frequency {l/tg ^ /o) to obtain better 
frequency resolution and reduce the statistical error. We 
then have 120 _(= T^s/ts) f^cq(/o) using Eq. Then, 
the estimator r2eq(/o) is 2.9j;Q24 x 10^^ from the expo- 
nential fitting of the r2oq(/o) distribution, rejecting the 
segments whose rioq(/o) are five times larger than the 
median calculated from all segments (see Fig|3]). The 
addition term is the statistical (fitting) error. 

The upper hmit on GWBs, i}^^{fo), can be described 
using a confidence level C defined by, 

r PinMm^^{fo))dn,,{fo), (4) 

where P(rics(/o)|^^^(/o)) is a conditional probability 
distribution and obeys the Gaussian distribution: 



Pi^csifa)\^gwifo)) oc exp 



(r!e.(/o)-r!gw(/o))' 



2n^ 



^gw(/o)V^ 

"(5) 

Here, ftcsifo) is the estimator of figw(/o) using the N{= 



o 4+0.26 . 
•^•^-0.29 



109) samples. As a result we found iiy^(/o) 
10" at C = 0.95. 

The biggest systematic error is the sensitivity from Al 
to 6. It arises from the uncertainty of the beam spot posi- 
tions in the end mirrors. From the size of the mirrors, we 
find that the maximum error is about 10%. For r2g^(/o), 
this error is about 20%, since ftcqifo) is proportional to 
^^(/o)- Considering the statistical error and the system- 
atic error, we finally obtain a conservative upper limit of 
rjUL ^ 4.3 X 1017 at C = 0.95. 
Discussion.- 

Comparison with the other experiments is shown in 
Fig. m Using the small-scaled TOBA, we successfully 
found the first nonintegrated upper limit on GWBs at 
the unexplored frequency that bridges the gap between 
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FIG. 5: Upper limits Red line is our new upper limit. 

Red dotted lines are our expected limits in next three phases. 
Current upper limits are also described as solid lines 
Dashed lines are the frequency-integrated upper limits [ll, l^- 



h ~ 10~^® Hz~^/^ at 0.1 Hz. Estimated upper limits 
at each phases are described in Fig. [5] with a one- year 
observation by a pair of two TOBAs. 

Conclusion- 

A TOBA has been proposed to search low-frequency 
GWs even with a ground-based configuration. We have 
developed a small-scaled TOBA and placed the first non- 
integrated upper limit on GWBs at 0.2 Hz. The new con- 
straint is r2gw(/) < 4.3 X 10^^ at 0.2 Hz with a bandwidth 
of 10 mHz. We experimentally demonstrated TOBA's 
capabilities and successfully opened the unexplored fre- 
quency band that the current GW antennas can not ac- 
cess. 

This work was supported by a Grant-in- Aid for JSPS 
Fellows, Young Scientists (A), GCOE for Phys. Sci. Fron- 
tier, MEXT, Japan. The authors thank I. Buder for care- 
fully reading the manuscript. 



the LIGO band (around 100 Hz) and the Cassini band 
(10~^ — 10"'' Hz). Considering integrated upper lim- 
its, the energy density r2gw(/) at 0.2 Hz is already con- 
strained by the BBN [l[ or the CMB measurements Q. 
However, astrophysical GWBs would be generated at 
much later times and, thus, would not be subject to the 
above limits. On the other hand, our limit can constrain 
such GWBs at 0.2 Hz. Therefore, our result complements 
other nonintegrated upper limits at different frequencies 
and integrated upper limits at 0.2 Hz. 

To buih the TOBA with the final configuration [13], 
we will have three phases. Applying well designed mag- 
netic shields to the small-scaled TOBA, the initial phase 
(phase-I) aims at the thermal noise limited noise level 
{h ~ 10^12 Hz-i/2 at 0.1 Hz) under lO"'^ Pa. Remember 
that the current thermal noise is limited by the residual 
gas noise. In this phase-I, we will move the small-scaled 
TOBA environmentally quieter site (Kamioka mine) and 
install the second test mass to obtain the common mode 
noise rejection for the seismic noise. The next phase 
(phase-II) is a middle-scaled TOBA using two 2 m scaled 
test masses. From the increase of the moment of inertia 
and the optimization of the mass-shape, force noises will 
be suppressed by a factor about 500. Target noise level is 
about 10-15 Hz-i/2 at 0.1 Hz using locked Fabry-Perot 
interferometers as the rotational sensors. After that, 
we will be in final phase (phase-Ill). A major change 
from the previous phase is 10 m scaled and cryogeni- 
cally cooled test masses. The cooling aims to suppress 
the thermal noises of the test masses and the suspen- 
sions. The TOBA with the final configuration will have 
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